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TWO-DIMENSIONAL DIFFUSION THEORY ANALYSIS OF REACTIVITY 
EFFECTS OF A FUEL-PLATE-REMOVAL EXPERIMENT 1 

By Edward R. Gotsky, James 1\ Cusick, and Donald Bocaut 


SUMMARY 

Two -dime nsiona l two-group diffusion calculations 
were performed on the NASA reactor simulator in 
order to evaluate the reactivity effects of f uel plates 
removed successively from the center experimental fuel 
element of a seven- by three-element core loading at 
the Oak Iiidge Bulk Shielding Facility. The reac- 
tivity calculations were performed by two methods: 
In the first , th e slowing-doum propert ies of the experi- 
mental fuel element were represented by its infinite 
media parameters; and f in the second, the finite size 
of the experimental fuel element was recognized, and 
the slowing-down properties of the surrounding core 
were attributed to this small region. The latter cal- 
culation method agreed very well with the experimental 
reactivity effects; the former method underestimated 
the experimental reactivity effects. 

INTRODUCTION 

Tlic NASA Lewis Research Center is construct- 
ing a high-flux research reactor with which the 
components of high power density reactors under 
development for various flight applications can be 
studied. The reactor employs an array of alumi- 
num fuel elements similar to those of the Materials 
Testing Reactor in Idaho. 

The geometric center of the active lattice of a 
reactor is an attractive location for experiments 
because of high flux and flux symmetry. The geo- 
metric center of the core, however, is also the most 
sensitive reactivity region for accidental composi- 
tional changes, and great care must be exorcised 
in the design of such in-pile experiments. Many 
reactivity calculations have been made for both 
fueled and unfueled experiments in a center test 
hole of the NASA research reactor. Since reac- 


tivity calculations arc subject to uncertainty, an 
experimental program with the Bulk Shielding 
Reactor at Oak Ridge National Laboratory was 
initiated in order to evaluate the methods of 
analysis. In these Oak Ridge experiments, the 
core configuration of the NASA research reactor 
was mocked up within the limits of excess reac- 
tivity and materials available to the Bulk Shield- 
ing Reactor. The loading was such that: the re- 
activity effects of voids within the core could he 
measured and the worth of fuel-element water 
passages and fuel-element plates could be deter- 
mined. The Bulk Shielding Reactor is fully dis- 
cussed in reference 1, and the experimental data 
for this particular loading as obtained at the Oak 
Ridge National Laboratory are presented in ref- 
erence 2 (pp. 8-10). 

The experiments in which the reactivity effects 
are measured are conveniently analyzed by group 
diffusion calculations. Within the limits of dif- 
fusion theory, two- and three-group analyses yield 
reasonable neutron flux distributions for thermal 
and epithermal reactors. F or hydrogen-moderated 
reactors, many modified group diffusion methods 
have been used to interpret criticality experiments 
(refs. 3, 4, and 5). 

The present report is concerned with the Bulk 
Shielding Reactor fuel-plate-removal experiments 
and the correlation of these data with diffusion 
theory calculations. These experiments have been 
analyzed by two-group two-dimensional diffusion 
calculations. The solutions have been obtained 
on a two-dimensional nuclear-reactor simulator 
that is based on the design of the simulator de- 
scribed in reference 6. 

The cooperation of the Oak Ridge National 
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Laboratory and the- staff of the Bulk Shielding 
Facility in performing the reactivity experiments 
is gratefully acknowledged. The authors are 
greatly indebted to Aliss Dorothy Hood, who 
supervised the calculations on the reactor 
simulator. 

REACTOR AND REACTIVITY EXPERIMENTS 

The Bulk Shielding Reactor is an assembly of 
fuel elements that may be arranged into various 
critical configurations. The fuel is highly enriched 
uranium contained in aluminum-clad fuel plates. 
A complete fuel element is made up of 18 fuel 
plates and contains a total of about 140 grams of 
uranium-235. 

REACTOR LOADING 

The reactor configuration used in the reactivity 
experiments is shown schematically in figure 1. 
The loading consisted of 21 fuel elements in a 
seven-by-three array with beryllium oxide (BeO) 
reflector pieces arranged in one row on the north 
and two rows on the south. The east and west 
faces of the core wore reflected by water ami per- 
mitted complete insertion or complete withdrawal 
of the two guillotine safety blades. Each safety 
blade consisted of a thin cadmium sheet between 
two aluminum plates approximately 12 by 24 


Number 2 guillotine 
safety blade--, 



Figure 1.— Reactor loading for reactivity experiments at 
Bulk Shielding Reactor. 


inches. The guillotine safety blades were guided 
and positioned, relative to the core, by grooved 
aluminum pieces and were supported by electro- 
magnets actuated by scram circuits. A regulating 
rod and a core safety rod were provided in the 
control-rod fuel elements in grid positions 23 and 
27, respectively. The regulating rod was cali- 
brated and used to measure the reactivity effects 
of core changes. 

The elements through which the conventional 
safety rod and the regulating rod move contained 
half the normal number of fuel plates and, there- 
fore, about 70 grams of uranium-235 each. The 
primary-reflector elements wore hot-pressed beryl- 
lium oxide blocks encased in watertight aluminum 
cans of the same outer dimensions as the fuel 
elements. The reactor was moderated and cooled 
by water that also served as the secondary reflector 
and reactor shield. The center grid position 25 
was occupied by the experimental fuel assembly. 

It should be noted that this core loading had 
vertical symmetry about the horizontal midplane 
and, in addition, east-west symmetry about a 
vertical plane across the narrow core dimension. 
(Although the top and bottom reflectors wine not 
identical, the reactor was symmetrical nuclear- 
wise because the top and bottom reflectors were 
predominantly water.) The vertical symmetry 
permits a typical horizontal slice through the 
reactor to be analyzed. Because of symmetry 
within the plane of this horizontal slice, only one- 
half oi this slice must be simulated for diffusion 
calculi lions. An additional calculation indicated 
that the effect of the partially inserted con' control 
rod dees not significantly alter the flux distribu- 
tions at the center grid position (25). 

REGULATING-ROD CALIBRATIONS 

In order to measure the reactivity effects of a 
wide range of compositional changes in the 
experimental fuel assemblies, it is necessary to 
have mi appreciable excess reactivity in the core 
loading. From safeguard considerations, core 
loadings at the Bulk Shielding Reactor are limited 
to excess reactivities of less than 2.5 percent. 
The regulating rod, which was used to measure 
the reactivity effects, must be worth less than 
1 percent. It was therefore calibrated by the 
method of distributed poisons for successive core 
configurations with both guillotine safeties out, 
one inserted, and finally both inserted. (Th 
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guillotine safeties were found to be worth about 
1 -percent reactivity each.) 

As the regulating rod in the 23 grid position 
was inserted, the safety rod in the 27 grid position 
was withdrawn to maintain criticality for the 
various guillotine configurations. This procedure 
was repeated for the clean core and for the core as 
uniformly poisoned as possible by 224 quarter- 
grain gold foils of 5-mil thickness. Seven gold 
foils were taped to each of 32 thin strips of Lueite. 
The foils were equally spaced on the strips and 
located in fuel-element water passages so as to 
poison the core vertically and laterally as uni- 
formly as possible. During these calibrations, a 
standard fuel-element assembly occupied the 
25 grid position. 

From elementary reactor theory, the change in 
reactivity resulting from the addition of the gold 
foils was taken as the fractional change in total 
macroscopic thermal absorption cross section of 
the core. Calibration of the regulating rod on 
this basis assumes that the gold foils act only as 
thermal absorbers. However, the strong reso- 
nance absorption of gold and the relatively large 
resonance flux in the core make epithermal ab- 
sorption by the gold significant. Therefore, the 
regulating-rod movement compensating for the 
gold foils is the result of thermal and epithermal 
absorption, and reactivity effect of the gold foils 
is only partially given by the fractional change in 
the total macroscopic thermal absorption cross 
section of the core. 

Although all the data presented herein and in 
reference 2 are based on calibrations neglecting 
resonance absorption in the gold foils, an estimate 
of the reactivity effects of resonance absorption 
may be made by using data presented in reference 
7. In reference 7, the reactivity worths of samples 
of gold of various thickness wore compared to the 
worth of samples of boron at various positions of 
the Reactivity Measurement. Facility, a low- 
power critical reactor located at the site of the 
Materials Testing Reactor at Idaho. The worth 
of boron samples in the fuel elements of this 
facility should represent the reactivity effects of 
a pure thermal absorber in a reactor like the Bulk 
Shielding Reactor reasonably well since the com- 
position and fuel loading are quite similar. Of 
course, the relatively larger worth of gold foils 
compared to boron for samples of the same total 
macroscopic thermal cross section, which is due 


to resonance absorption at the surface of the gold 
samples, will decrease with thicker gold foils 
because of the self-shielding of the resonance 
neutron flux. Indeed, the data of reference 7 
indicate the ratio of reactivity of gold to boron of 
comparable total macroscopic thermal cross sec- 
tion to vary from 1.6 for gold-foil thickness of 1 
mil to 1.3 for gold-foil thickness of 5 mils. The 
effect of using a correction factor of 1.3 in the 
experimental reactivity effects is illustrated later. 
Because there is some epithermal absorption in 
boron, the gold-foil correction factor is probably 
somewhat greater than 1.3. 

REMOVABLE-FUEL-PLATE ASSEMBLY 

In the reactivity experiments discussed in this 
report, grid position 25 was occupied by a standard 
fuel element in which only 6 of the 18 curved fuel 
plates were brazed into the grooved side plates. 
A photograph of this removable-fuol-plate assem- 
bly is shown in figure 2. The fuel plates are 
removable with the exception of the two end 
plates and the four plates blocked by the fuel- 
element handle. The plates are 60 mils thick, and 
the water passages are 117 mils thick; the fuel 
element is 24 inches long and 3 inches wide. 

These fuel plates were removed in two sequences. 
In sequence I, the interactions were maximized by 
successively removing adjacent fuel plates and 
thereby forming an enlarging water region near the 
center of the cone. In sequence 11, the inter- 
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actions of each plate removed were initially 
minimized by removing plates farthest apart from 
each other. These sequences are illustrated in 
figure 3, in which letters A to L represent the 
fuel plates removed in each step. 

Standard fuel element 


12 Removable 

fuel plates- ABODE FGHIJKL 

Figure 3. — Sequences of removing fuel plates in standard 
fuel element of Bulk Shielding Reactor. 

EXPERIMENTAL RESULTS 

The reactivity per incremental fuel plate re- 
moved in the maximum interaction (sequence 1) 
and in the minimum interaction (sequence II) is 
presented in figure 4. The net effect of the re- 
moval of fuel plates and the consequent enlarge- 
ment of water regions decreased the reactivity in 
all cases. The reactivity for eacli fuel plate re- 
moved was greatly affected by the sequence of 
removal. In sequence I, for example, the worth 
of adjacent fuel plates increased rapidly as the 
region devoid of fuel enlarged. On the other 
hand, in sequence II, individual fuel plates from 
opposite sides of the element had equal worth, 
while the worth of fuel plates from intermediate 
positions in the element increased slowly at first. 

In these rcmovable-fuel-plate sequences, avail- 
able core excess reactivity permitted 11 fuel plates 
to be removed. In the maximum interaction, 
removal of the twelfth plate (number D in fig. 
4(a)) shut down the reactor when the plate was 
45 percent withdrawn from the core. At this 
point, the regulating rod was fully withdrawn 
from the core; this indicated that the withdrawn 
portion of the twelfth fuel plate was worth 
0.16 AKjK percent. (All symbols are defined in 
appendix A.) The accumulated reactivity of the 
first 11 fuel plates was 2.36 percent. The com- 
plete twelfth fuel plate is estimated to be worth 
at least 0.16/0.45=0.36 AKjK percent (shown in 
fig. 4(a) as a tailed datum point). In the mini- 
mum interaction, the worth of the twelfth fuel 
plate is readily estimated by subtracting the 
cumulative reactivity of sequence II, step 11, 
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Fig i IE 4. Reactivity for each fuel plate removed. 


from t ie cumulative reactivity of sequence I, 
step 12 A reactivity of 0.48 AKjK percent was 
obtained, as shown in figure 4(b) (tailed datum 
point). 

The 'umulative reactivity of these sequences 
is prose ited in figure 5 as a function of the number 
of fuel plates removed (each fuel plate contains 
7.78 g <>;* uranium-235). The minimum-interaction 
sequent e was used as a basis for the reactivity 
calculations reported herein and for the comparison 
betweei the analytical group diffusion results and 
the experimental data. 


GROUP DIFFUSION CALCULATIONS ON REACTOR 
SIMULATOR 

SIMULATOR 

The NASA two-dimensional reactor simulator is 
an anal )g computer for obtaining solutions to the 
neutron -group diffusion equat ions. Reactivity 
problems in which spatial symmetry exists in only 
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Figure 5. — Cumulative reactivity of fuel plates in element 
at grid position 25. 

one dimension may be solved. The simulator is 
a direct extension of the one-dimensional device in 
use at the NASA Lewis Research Center since 
1951 (ref. 6). 

The two-dimensional reactor simulator is shown 
in figure 6. The fast- and thermal-neutron-group 
resistor-network boards, in which the space -point 
resistors for a particular problem are prefabricated 
in aluminum plug-in units, may be seen on the 
right. The source-input panels are shown on the 
left. The problem is iterated from group to group 
by varying the source input until a converged 
flux distribution and value of the effective multi- 
plication factor K ef f are obtained. In the two- 
group formulation used, epithermal fission and 
absorption are part of the fast-group processes. 



Figure (>. — Two-dimensional reactor simulator. 


SIMULATOR REPRESENTATION OF LOADING 

The core loading shown previously had vertical 
symmetry about the horizontal midplane and, in 
addition, symmetry about a vertical plane across 
the narrow core dimension. The vertical sym- 
metry permits a typical horizontal slice through 
the reactor to he analyzed. Because of symmetry 
within the plane of this horizontal slice, it is neces- 
sary to simulate only half of this slice for the 
diffusion calculations. Representation of this 
geometry on the reactor simulator, therefore, 
takes the form shown in figure 7, where the core 
fuel elements, the primary beryllium oxide reflec- 
tors, and the secondary water reflector are viewed 
from above. The fuel element in the central 
core position, from which fuel plates were succes- 
sively removed, is shown in the inset. There is 
a net point for each area shown. Each coarse 
net point has the dimensions of a fuel element, 
about 3 inches square. The medium-sized net 
points are one-fourth of the coarse grid size. The 
area at the center of the core represents one-half 
of the experimental fuel element from which fuel 
plates were individually removed. In order to 
obtain a detailed description of the fluxes in this 
region, finely sized net points that represent 
areas about a centimeter on a side were used. In 
all, the calculation involved 69 coarse net points, 
34 medium net points, and 32 fine net points in 
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Removable-fuel-plate assembly at center of core 

Figure 7- -ltcactor geometry for two-dimensional 
simulator solutions. 
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each group of the two-group solution. The 
nuclear properties used for each space point are 
for the homogenized contents of the experimental 
region represented. 

TWO-GROUP FORMULATION 

The two-group multiregion flux equations and 
a brief description of the methods of evaluating 
group parameters are given in appendix B. The 
two-group constants used in these calculations 
for the unperturbed core are shown in table I. 
These parameter values are similar to those used 
to predict reactivity effects in the early Materials 
Testing Reactor (Idaho) loadings, where fuel 
elements similar to those in the Bulk Shielding 
Reactor were used. In the two-group formulation 
used here, epithermal fission and absorption are 
retained as part of the fast-group processes; this 
requires the use of the fast multiplication constant 
K f and the nonabsorption probability B th . It is 
shown in reference 8 that a two-group formulation 
including epithermal absorption and fission as used 
herein predicts reactivities that agree with a 
formulation which considers in detail the energy- 
dependent nuclear processes in the epithermal 
region. Values of L 2 and L 2 th are derived from 
experimental values for metal-water mixtures; and 
other nont hernial values are calculated from 
si owing-down distributions in infinite media re- 
gions of the same composition. Parameter values 
for regions in the center fuel element with fuel 
plates removed were estimated in this manner also. 

In all eases, the neutron leakage in the direction 
normal to the plane of the two-dimensional solu- 
tion is assumed proportional to the local flux with 
the vertical geometric buckling B\ as the constant 
of proportionality. The B 2 value of 0.001(122 is 
based on an actual core height of 02.0 centimeters 
plus reflector savings of 10.0 centimeters due to 


TABLE I.— 1 TWO-GROUP DIFFUSION 
PARAMETERS 

(Vertical buckling in all cases, R*=^O.OOI622.J 
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the aluminum-water regions above and below the 
core. The problem of keeping this vertical buck- 
ling constant as more and more fuel plates are 
removed accounts for part of the small discrep- 
ancy between experiment and calculations. 

The question immediately arises as to the va- 
lidity of diffusion theory in regions that are geo- 
metric illy small, such as regions where a few 
fuel plates have been removed. For example, 
removal of two adjacent fuel plates results in a 
water passage that is about 1 by 7 centimeters 
in cross section and runs the full core height. 
The slowing-down properties of such a small region 
can In rdly be characterized by the slowing-down 
length in water, which is about 5.0 centimeters. 
It is reasonable to assume that the slowing-down 
properties of the surrounding core may be more 
significant. For this reason, the diffusion calcu- 
lations for each core configuration were performed 
twice: first, using the homogeneous diffusion prop- 
erties for each region regardless of size (referred to 
hereafter as “local region slowing-down” cases); 
and second, using the values of L } and p th for the 
surrounding core to characterize the slowing down 
with absorption in the regions where fuel plates 
were removed (referred to as “surrounding region 
slowing-down” cases). In the surrounding region 
slowing-down cases, two-group constants other 
than and p th were based on the local region 
proper t ies. 

EXPERIMENTAL CONFIGURATIONS CALCULATED 

( 'all illations of the reactivity effects resulting 
from fuel-plate removal were performed for the 
referoi ce core with no fuel plates removed and 
for tin minimum interaction sequence in which 6, 
10, ai d the equivalent of 13.5 fuel plates were 
removed from the center fuel element. The fine- 
grid-point arrangement used in each case dupli- 
cated 'jompositionally the experimental geometry 
within the limits of the 32 space points available. 
An additional case, in which the range of the 
experimental data was extended by removing all 
18 plates of the central fuel element, was calcu- 
lated i i order to estimate the worth of one complete 
fuel element at the center of the core. 

CALCULATED RESULTS 

UNPERTURBED CORE 

The representative two-dimensional flux dis- 
trihut ons are shown in figure 8 for the unper- 
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turbed core with no fuel plates removed from the 
center fuel element. Figure 8(a) shows the inte- 
grated fast-llux distribution in the core and re- 
flectors. The magnitude of the fluxes is shown 


relative to a spatial average core thermal-neutron 
flux of unity. Figure 8(b) shows the thermal- 
neutron Uux distribution with characteristic* peak- 
ing in the reflectors and falloff toward tin* cornel's 
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of tlip slab. There is an over-all asymmetry due 
to the thicker primary reflector on one face of the 
core. 

In order to observe the effects of fuel-plate 
removal from the central fuel element, the flux 
distributions across the narrow core dimension 
where maximum perturbations occur are noted. 
Reactivity effects were evaluated by integrating 
the two-dimensional flux distributions obtained 
over the core volume. 

On an absolute basis, the value of K tff for this 
Bulk Shielding Reactor core loading was observed 
experimentally to be 1.025; that is, the loading 
had 2. 5-percent- excess reactivity. The value of 
K €fI estimated from two-dimensional diffusion cal- 
culations is 1.032, an estimated excess reactivity 
of 3.2 percent. 


LOCAL REGION SLO WING-DO WN MODEL 

Neutron flux distributions for several local 
region slowing-down cases are shown in figure 9. 
The flux distributions across the narrow core 
dimension art* shown for three cases: (1) the un- 
perturbed con*, (2) 10 fuel plates removed, and 
(3) all 18 fuel plates in the center fuel element 
removed. The fluxes are relative to a spatial 
averag * thermal flux of unity in (lit* core. Shown 
in figu v 9 are data for the core with the* central 
fuel element, the beryllium oxide side reflectors of 
different thicknesses, and the secondary water 
reflectors. Removal of the fuel plates and, there- 
fore, fast-neutron sources, with consequent re- 
placement. by water, lias reduced the fast flux 
tp f for it large region in the vicinity of the central 
fuel element, by providing a large internal slowing- 
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down region. The thermal fluxes <p th arc conse- 
quently limit up in this region. The overflow of 
mill Irons diffusing into surrounding fuel elements 
is felt, for about two thermal-diffusion lengths into 
the core, after which the thermal fluxes are 
unperturbed. The reactivity worth of remaining 
center fuel plates and of plates in adjoining ele- 
ments is increased as a result. Tin* net reactivity 
effect- is the sum of the —A K due to fuel-plate 
removal and the + AK due to adjoining fuel 
plates operating in higher thermal fluxes. 

SURROUNDING REGION SIX) WING-DOWN MODEL 

Figure 10 illustrates the two-group neutron 
fluxes for 10 fuel plates removed, as calculated by 
assuming the values of Lj and j) th for the missing- 
fuel-plate region to be the same as those of the 
surrounding core. These flux distributions are 


compared with the fluxes for the unperturbed 
core. The slo wing-down parameter L 2 f for the 
core is larger than the value of TJ f for the missing- 
fuel-plate region, with the result that fewer 
neutrons slow down in this central region. Con- 
sequently, the thermal-neutron flux at the center 
region for this surrounding region slowing-down 
case is lower than the thermal-neutron flux for the 
corresponding local region slowing-down case. 
The remaining fuel plates in the center test ele- 
ment and in the adjoining elements, being in 
regions of lower thermal flux for the surrounding 
region slowing-down case, are worth less than 
the fuel plates for the local region slowing-down 
case where thermal flux is higher. Inasmuch as 
the net reactivity effect is the sum of the —AK 
due to fuel-plate removal and the + A/v due to 
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adjoining fuel plates operating in higher relative 
thermal fluxes, quite different effects may he 
expected from the two diffusion calculations. 

DISCUSSION 

The net reactivity effects of fuel-plate removal 
from the center fuel element, as measured at the 
Bulk Shielding Reactor and as presently calcu- 
lated, are compared in figure 11. Two experi- 
mental broken curves are shown; the one with 
the circular data points presents the data from 
the calibration curve in which the resonance 
absorption in the calibrating gold foils was 
neglected. The other broken curve has been 
corrected by a factor of 1 .3 as discussed previously. 

The local region slowing-down calculations 
consistently underestimate the observed reac- 
tivity effects. On the other hand, the surrounding 
region slowing-down calculations predict reac- 
tivity effects for fuel-plate removal up to about 
eight plates quite well. 

Both calculations underestimate an extrapola- 
tion of the experimental data corresponding to 
removal of the remaining fuel plates in the center 
fuel element. This may indicate that a third- 
dimensional effect is becoming important. In all 
the present calculations, neutron leakage in the 
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Fict’itK 11. Comparison of experimental and calculated 
reactivity effects of fuel-plate removal. Bulk Shielding 
Reactor experiments, minimum-interaction sequence. 


vertica (z) direction, that is, normal to the plane 
of solution, was taken to be proportional to the 
local flux ^ with a fixed value of the vertical 
geometric buckling B\ as the constant of propor- 
tionality. This assumes a cosine flux distribu- 
tion in the ^-direction above and below the plane 
of solution. Furthermore, this assumption im- 
plies knit the flux in any horizontal plane is 
dh ■ectly proportional to the flux in the horizontal 
midpla le. 

As t io missing-fuel-plate region at the center 
of the core enlarges, however, an increasingly 
important diffusion region with reflector proper- 
ties is provided; this gives rise to a nonuniform 
flux distribution above 1 and below flu* plane of 
solutio n This internal water region acts essen- 
tially is a fast control rod by collecting fast 
neutrons over a relatively large range of several 
fast diffusion lengths in the core, by moderating, 
and then by diffusing them hack into the core 
as thermal neutrons over a shorter range of several 
thermal diffusion lengths. 

To illustrate 1 his third-dimensional effect, con- 
sider figure 12, in which portions of the center 
fuel element are removed from the top and bot- 
tom of the eon* and replaced with water. This 



Fun kk 12. Illustration of asymmetry due to water 
passage in core center. 
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internal reflector (water) gathers in fast neutrons 
from all directions over a range of about one fast- 
transport mean free path in the core (4 cm) and 
moves the neutrons away from the center of the 
core. These neutrons are thermalyzcd in the 
water region and diffuse into the core over a 
range of about one thermal-transport mean free 
path (0.8 cm). The net result is a transfer of 
fast neutrons away from the center of the core 
to regions of less importance. As the internal- 
reflector region increases in height, the migration 
of fast neutrons away from the core center also 
increases. 

The two-dimensional diffusion calculations can- 
not account, for this fast-neutron migration away 
from the center of the core, and therefore higher 
thermal fluxes are indicated in the central core 
region than are present in the experiment. The 
reactivity worth of the fuel elements surrounding 
the central fuel element is consequently overesti- 
mated, and the worth of the central fuel element 
is underestimated. 


An approximation of the diffusion in the in- 
ternal water region can be made by decreasing 
the effective average height of the core; this re- 
sults in an increase in B\ and a consequent increase 
in the vertical leakage term D f B 2 z . To illustrate 
this, a reduction in effective core height of 1 
centimeter results in a decrease in reactivity of 0.3 
percent A K for this core loading. (The reduction 
in reactivity is linear for about 6 cm in effective 
core height.) 

CONCLUSION 

The three-dimensional reactivity effects meas- 
ured in the Bulk Shielding Reactor core-loading 
experiments are approximated by a two- 
dimensional calculation. Over the range of 
experimental data, the two-dimensional diffusion 
calculations predict reactivity effects that agree 
with the experimental data reasonably' well. 

Lewis Research Center 

National Aeronautic* and Space Administration 
Cleveland, Ohio, October 18, 1957 



APPENDIX A 


SYMBOLS 


m 

buckling in direction normal to plane of 

5 1 

top reflector savings 


solution 

b 1 

bottom reflector savings 

D 

diffusion coefficient 

V 

" A 

macroscopic absorption cross section 

h 

height of active lattice 

v 

^ F 

macroscopic fission cross section 

K 

multiplication constant 

V 

~Sl 

macroscopic slowing-down cross section 

K ef r 

effective multiplication factor for reactor 

V 

number of neutrons emitted per fission 

Kj 

mean-square slowing-down length for fis- 
sion neutrons 

<p 

neutron flux 

L,\ 

mean-square diffusion length for thermal 
neutrons 


SUBSCRIPTS 

pth 

fraction of neutrons escaping capture in 

f 

fast neutron group 


slowing down 

th 

thermal neutron group 
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APPENDIX B 


DIFFUSION EQUATIONS AND GROUP CONSTANTS USED IN REACTIVITY CALCULATIONS 


EQUATIONS 

Two-group multiregion two-dimensional diffu- 
sion calculations are based on the following flux 
equations that include effects of fast absorption 
and fission: 


(*+*') 

\dx 2 ^dy 2 J 




+ J'Sjf-./wf th<Pth = 0 

(Bl) 

( d 2 d 2 

) <Ptn— (Sa, 



(B2) 


where x and y are the coordinates in the plane of 
solution. In equations (Bl) and (B2), effective 
absorbers are introduced to account for neutron 
leakage in the direction normal to the plane of 
solution; these are D f B\ and D {h B 2 for the epi- 
thermal and thermal groups, respectively. The 
quantity B\ is the geometric buckling given 
in terms of the equivalent bare-core dimensions 
in the direction normal to the plane of solution. 
For the present reactor, 

B 2 = 

2 (h+5 1 +5 2 )* 


The multiplication constants representing the 
number of neutrons born per neutron absorbed 
in each group are 

Kf=y (B6) 

^A,f 

K th =v^ (B7) 

With these definitions, two-group equations 
(Bl) and (B2) may be rewritten as: 

+K, (1 -p th ) % <Pf+K,n yv *»=0 (B8) 

-I'f 

D " (dx 2+ dy 2 ) Vtk+P,h V, 

(B9) 

Equations (B8) and (B9) will apply for passive 
regions if K f —K ih — 0. At the interfaces between 
regions, the fluxes and neutron currents are con- 
tinuous. At the extrapolated boundaries of the 
secondary reflector, the fluxes are zero. 


Inasmuch as the macroscopic cross sections 
are the effective averages in each energy group, 
the mean-square slowing-down length L 2 f and the 
mean-square thermal diffusion length JJ th may be 
defined as 


T 2 _ Df 

1 ^a, 




D tl 


(B3) 

(B4) 


The fraction of neutrons escaping capture in 
slowing down also may be expressed in terms 
of these effective values as 


Pth 


^Sl,f 


(B5) 


GROUP PARAMETERS 

The constants required for solution of the two- 
group equations (B8) and (B9) in each region are: 

(1) Fast: L 2 n D h K f , p th 

(2) Thermal: L* th , D th , K th 

The validity of a group representation depends 
greatly on the correctness of the group constants 
used to represent the various competing nuclear 
processes. Procedures used to evaluate group 
constants for highly enriched uranium hydro- 
geneous media are described in references 5 and 
8 and are outlined in the following section. 

Fast constants.— For neutrons slowing down in 
hydrogeneous media, the bulk of the contribution 
to L/ is made at high neutron energies, that is, 
at energies above the epithermal region for which 
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the inverse velocity variation of microscopic 
absorption cross section has reduced absorption 
to a negligible value. For highly enriched ura- 
nium systems, most of the absorption occurs in 
the epithermal region in the energy range below 
about 1000 electron volts. Hence, fast-group 
diffusion, leakage, and slowing-down processes 
may be effectively separated from epithermal 
absorption processes. If the Fermi age theory is 
taken to apply in the epithermal region, the resid- 
ual slowing-down, diffusion, and absorption proc- 
esses in this region may be treated in detail. To 
accomplish this, neutrons should be divided into 
three separate groups: fast, epithermal, and 

thermal, as in references f> and 8. All neutron 
absorption, and therefore production, is restricted 
to the epithermal and thermal regions. Group- 
diffusion slowing down in the fast region is then 
characterized by values of L) estimated from cor- 
relations of the few experimental metal-water- 
mixture measurements (ref. 5). Slowing down in 
the epithermal region is calculated by the Fermi 
age theory, which provides the variation of neu- 
tron flux with energy to permit evaluation of the 
effective values of X A , and#*/,. These average 
values are obtained by weighting local values by 
the neutron flux in the epithermal region. 

In the present two-group calculations, epither- 
mal fissions and absorptions have been retained 
in the fast group, and the age-theory neutron 
slowing-down distribution in infinite media has 
been used to evaluate the fast-group constants. 
The fast constants D fl K fy and j p lh are, therefore, 
obtained by weighting local values according to 
the energy distribution of neutron flux in an 
infinite medium of the same composition as the 
pertinent reactor region and as indicated by age 
theory. The fission spectrum is included in de- 
termining the energy distribution of neutron flux. 


It is shown in reference 8 that a two-group 
formulation including epithermal absorption and 
fission predicts reactivities that agree with a 
formulation which considers in detail the energy- 
dependent nuclear processes in the epithermal 
region 

Thermal constants. — The thermal parameters 
are calculated from equations (B4) and (B7). 
Values of th and 'Z Ft th are macroscopic cross 
sections obtained by averaging local values over 
the Maxwellian distribution of thermal neutron 
flux. The effects of thermal spectrum hardening 
due to preferential absorption in the lower energy 
portions of the thermal distribution have been 
neglected. 

The diffusion coefficient D th is given by the sum 
of tin contributions of each constituent atom, 
the hydrogen atom usually being the largest con- 
tributor. The effects of the chemical binding of 
the hydrogen atom to the water molecule, how- 
ever, significantly alter thermal scattering prop- 
erties of hydrogen because of the increasing effec- 
tive mass of the hydrogen atom for neutrons in the 
lower energy portions of the thermal distribution. 

The effective value of D lfi for any media is also 
obtained by averaging local values over the 
Maxwellian distribution of neutron flux, provided 
the variation of effective mass of the hydrogen 
atom with local energy is known. In the absence 
of these specific data, a method attributable to 
Radkowsky (described in ref. 9), which checks 
experimentally determined thermal diffusion prop- 
erties of water for a range of temperature, was 
used. A value of 31.4 barns was obtained for the 
effective microscopic transport cross section for 
hydrogen. This value, which corresponds to an 
experimental value of thermal diffusion length 
of 2.89 centimeters in water at room temperature, 
has been used in the present calculations. 
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